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Abstract: An improved equivalent simulation model for a CMOS -integrated Hall plate is 
described in this paper. Compared with existing models, this model covers voltage 
dependent non-linear effects, geometrical effects, temperature effects and packaging stress 
influences, and only includes a small number of physical and technological parameters. In 
addition, the structure of this model is relatively simple, consisting of a passive network 
with eight non-linear resistances, four current-controlled voltage sources and four parasitic 
capacitances. The model has been written in Verilog-A hardware description language and 
it performed successfully in a Cadence Spectre simulator. The model's simulation results 
are in good agreement with the classic experimental results reported in the literature. 
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1. Introduction 

Presently, CMOS integrated Hall magnetic sensors are widely used in many practical fields. 
Besides directly measuring the value of magnetic field, they are usually used to indirectly measure 
position, distance, speed, rotational angle or an electric current [1,2]. For instance, they can act as an 
automotive vehicle speed sensor, a replacement for mechanical switches, a brushless control for DC 
motors, and so on. Unfortunately, CMOS integrated Hall sensors have traditionally suffered from a lot 
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of non-idealities, such as low sensitivity, large offset, temperature drifts, non-linearity and packaging 
stress influence etc., which severely deteriorates their performance [3]. As a consequence, CMOS 
integrated Hall devices must depend on the processing circuit for offset and noise cancellation, 
temperature compensation and non-linearity correction. In order to facilitate the simulation analysis of 
electrical circuit with integrated Hall devices, it is necessary to extract a precise simulation model to 
take into account important physical effects and technological influences. Furthermore, the extracted 
model should be simple and conveniently implemented in standard SPICE- like EDA tools. 

Several compact simulation models of Hall elements have been put forward. Previously reported 
4-resistance Wheatstone bridge models don't fully take into account correlative physical and geometrical 
effects such as non-linear conductivity, junction effect, temperature drift, frequency-response, noise 
behavior and device shape-dependent sensitivity [4,5]. Later, Dimitropoulos et al. proposed a completely 
scalable lumped-circuit model to analyze all those effects, except for the influence of packaging 
stress [6]. The basic component for the lumped-circuit model consists of JFETs and current-controlled 
current sources. The number of these components can be freely increased to achieve the required 
accuracy at the expense of computation efficiency. However, this macro model needs an accurate 
JEFTs device model which normally cannot be provided by the standard CMOS technology. Recently, 
Madec et al. developed a compact model of a cross-shaped horizontal integrated Hall sensor [7]. It 
uses six sub-components to accurately model the non-linear resistance, allowing for the influence of 
space charge region modulation due to sensor bias. Unfortunately, it cannot consider sensitivity drifts, 
temperature drifts and influence of mechanical stress on offset. Besides, the resistance computation of 
the model has to be fed by empirical parameters through FEM (finite element method) simulation. 

In this paper, an accurate 8-resistance simulation model for a cross-shaped CMOS-integrated Hall 
plate is developed. To be conveniently used by circuit designers, this model is improved by replacing the 
JFETs with passive non-linear resistances and depletion capacitances. It takes into account voltage 
dependent non-linear effects, geometrical effects, temperature effects, and packaging stress influence, etc. 
Since we mainly deal with the magnetic sensors operating in a weak magnetic field in this work, two 
additional strong magnetic field related effects, namely magneto-resistance and carriers scattering, are not 
included in this model. The model has been written in Verilog-A hardware description language and was 
successfully tested in a Cadence Spectre simulator. This paper is organized as follows: in Section 2, we 
introduce the structure of the model and analyze the important physical effects of the Hall device with the 
basic equations. Furthermore, the detailed computation of device parameters in this model is presented. 
In Section 3, the simulation results of the model are compared with the classic experimental results 
reported in the literature. Section 4 summarizes this paper with some ideas for future work. 

2. The Improved Compact Model 

To be compatible with the spinning current techniques for reducing Hall offset [3], 90° symmetry 
Hall plates with square or cross-shaped structures are usually recommended. We can obtain the 
Z-matrix for the 90° symmetry Hall plates with four contacts illustrated in Figure 1 as follows [8]: 
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Figure 1. Diagram of measuring Z-matrix for a 90° symmetry Hall plate. 

1/2 




If the fourth contact is applied to the reference ground, the Z-matrix of the 90° symmetry Hall plate 
is only decided by three parameters Z x x , Z 12 , Z 13 . If the input current / t is applied to the first contact, 
the three measuring potentials shown in Figure 1 have the following relation: U l -U 2 =U 3 , and then 
we can obtain Z n — Z 12 = Z 13 . As a result, 90° symmetry Hall plates require at least two types of 
resistances to model their electrical properties. Thus, an 8-resistance model topology for the 90° 
symmetry Hall plate is suggested, which is illustrated in Figure 2. Its Z-matrix at zero-magnetic field is 
expressed by: 
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Figure 2. An equivalent model topology for the 90°symmetry Hall plate. 




However, in the conventional 4-resistance Wheatstone bridge model [4,5], the diagonal 
resistances R D are often neglected. The value of the resistance between two adjacent contacts is not 
accurate because current lines linking the two adjacent contacts do not flow across the center of the 
device. But this problem can be solved well in the 8-resistance model so that the simulation accuracy 
can be improved. 

2.1. The Structure of the Model 

The 90° symmetry cross-shaped Hall plate (see Figure 3) is most widely used because of its high 
sensitivity and immunity to alignment tolerances resulting from the fabrication process. Its fabrication 
technology is fully compatible with the standard CMOS process. As shown in Figure 3(a), the active 
area of the cross-shaped Hall plate is usually realized by a weakly doped N-well diffusion region. The 
N-well is isolated from the P-type substrate by the reverse-biased well/substrate p-n junction. A 
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shallow heavily doped top P+ layer often covers the surface of active area to decrease the flicker noise 
and the surface losses. It is normally formed to create the source and drain regions of PMOS 
transistors [9]. In addition, four contact regions are highly N+ doped to reduce the contact resistances in 
the source and drain formation processing step for NMOS transistors. Based on the basic model topology 
shown in Figure 2, a new simplified 8-resistance model for the CMOS integrated cross-shaped Hall plate 
is developed, as shown in Figure 4. Compared with Dimitropoulos' model, the simulation model is 
improved by replacing the JFETs by N-well body resistances and depletion capacitances. The four 
depletion capacitances are added into the model to simulate the transient behavior of the Hall plate. 
Besides, there are four controlled voltage sources V HI2 to model the Hall voltage. Each Hall voltage 
source V H/2 is controlled by the electrical current flowing through the nearer contact. 



Figure 3. Cross-shaped Hall plate fabricated in standard CMOS technology, (a) Top view; 
(b) View with the cross-section along B and D contacts. 
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In order to determine the resistance values of R H and R D in our model, a new and simple 
computation method is proposed in contrast to the FEM simulation method [7]. It is well known that it 
is best to measure the N-well sheet resistance 7? s according to the Van-der-Pauw method. Although the 
Van-der-Pauw method requires that the contacts of Hall device be point-like, it has been reported that 
cross-shaped Hall plate with a finger length to width ratio larger than 1:1 can give an accurate R s value 
with an error of less than 0.1% [10]. Usually, the required ratio of finger length to width can be 
fulfilled for achieving a high current related sensitivity. Thus, in the case of symmetric Hall plates, the 
sheet resistance can be determined by measuring the resistance value R AB CD in term of 
Van-der-Pauw method [11]: 
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where R m ^ = V CD //^presents the voltage difference between contacts C and D dividing the current 
flowing from contact A to contact B. The contacts of A, B, C and D are shown in Figure 3(a). 

Figure 4. A simplified model for the CMOS integrated cross- shaped Hall plate. 
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On the other hand, according to the structure of the model illustrated in Figure 4, the resistance 
r ab,cd i s calculated by: 
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The internal resistance between two diagonal contacts is given by: 
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Here, (2IJW+2/3) is the effective square number of the N-well resistance. L and W are the finger 
length and finger width of cross-shaped Hall plate, respectively. The center square number is 
approximately reduced to 2/3 as the two fingers for sensing Hall signal are placed in parallel. 
Considering the Equations (3), (4) and (5), finally we can obtain: 
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The N-well sheet resistance R s is calculated by: 
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Here, N D NW is the N-well doping concentration, t eff is the effective depth of Hall plate. As shown 
in Figure 3(b), it is equal to: 

teff = ^NW ~ tp+ ~ W NW,SUB ~ W NW ,P+ (9) 

where t NW is the depth of N-well implantation region, is the thickness of the top P+ layer, w m SUB 
is the bottom depletion region situated in between N-well and P-type substrate, w m p+ is the upper 
depletion region situated in between N-well and top P+ layer. 

Note that there are two main parasitic capacitances distributed across the Hall device body: (1) the 
reverse-biased upper depletion capacitance between the top P+ layer and N-well; (2) the 
reverse-biased bottom depletion capacitance between the N-well and p-type substrate. Usually the top 
P+ layer and P-type substrate are applied to ground together, thus they are connected in parallel 
physically. Unfortunately, the parasitic capacitances may limit the switching frequency for the 
spinning current offset reduction method. In order to simulate the complete ac behavior of the Hall 
plate, these parasitic depletion capacitances should be included in the model. Assuming the one-sided 
abrupt junctions, each depletion capacitance per unit area is calculated by following Equation (5): 
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Here, V hi is the built in potential of PN junction, N D NW is the doping of N-well, and N A is the 
doping of top P+ layer or P-type substrate. 

2.2. Hall Voltage and Magnetic Sensitivities 

When a magnetic field B is orthogonally applied on a device plane and two diagonal contacts are 
biased with a current /, the Hall effect takes place. Then the hall voltage V H appears on two additional 
contacts, it is equal to: 

V H =S I IB (11) 

where S { is the current related sensitivity. It depends on device geometry (geometrical correction 
factor G and Hall plate effective thick t eff ) and technology parameters (N-well doping N D NW and Hall 
factor r H ) [12]. It is defined as: 

H ly D,NW l eff 

When Hall plate is biased with a voltage source V, the Hall voltage is expressed by voltage related 
sensitivity S v : 

V H = S V VB (13) 

where S v = jU H G/ N square . ju H = r H /u n , it is the Hall mobility, //„ is the electron mobility, and N squre is 
the equivalent square number of N-well diffusion resistance between two diagonal contacts. 

The impact of the Hall devices geometry on Hall voltage is modeled by a geometrical correction 
factor. For a cross-shaped Hall plate, it can be calculated by using a conformal mapping [12]: 
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where 0 n = tan l (/u H B) , it is defined as the Hall angle. If W/2L < 0.39 , G has an accuracy better than 
0.5% [12]. 

In our model illustrated in Figure 4, each Hall voltage V H/2 is modeled by using the 
current-controlled voltage sources with the following equation: 

y H ,2 = ^S I I(n 1 ,n 2 )B (15) 

with /(np n 2 ) being current flowing between the contacts n l and n 2 . 

2.3. Voltage Dependent Non-Linear Effect 

It is well known that the thickness of depletion region is obviously changed by the reverse biased 
PN junction. Therefore, both sheet resistance and magnetic sensitivity suffer from a strong voltage 
non-linearity dependence. Since the doping concentration of the P+ top layer is obviously higher than 
that of the P-type substrate, the thickness variation of the upper depletion region modulated by reverse 
biased voltage can be approximately ignored. Using Equation (8), which is extended by the voltage 
dependent t eff , and a Taylor expansion results up to second order are given by [5]: 
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where R s (0V) is the zero-biased N-well sheet resistance, bbr^ and bbr 2 are the first and second 
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voltage dependency of resistance coefficients, respectively. k x = — — , and 

1 (N ASUB + N DNW )N DNW 

Cff = f Nw ~tp+ ~ w n,p+(QV) ■ Here, w N ^+(0^) is the upper depletion region with zero-bias. 
With the same calculation method, the current related sensitivity is modeled by: 

S^t/pJ = S^ + BBS x U pn +BBS 2 U 2 pn ) (17) 

where BBS l and BBS 2 denote the first and second voltage dependent coefficients of sensitivity, 
respectively. 

2.4. Temperature Effect 

We know that the temperature drift has serious effects on the equivalent N-well resistance, 
sensitivities and offset of the Hall device. The temperature behavior of N-well sheet resistance can be 
well approximated by the second order polynomial: 
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R s (U pn ,T) = R s (U pn ,30OK)-[l + R Tcr (T-30OK) + R TC2 -(T-300K) 2 ] (18) 

where R TCl and 7? rc2 are temperature coefficients of N-well resistance. These parameters can be directly 

obtained from foundry technological files. Rs(U pn ,300K) is the sheet resistance at the room temperature. 

Since the thermal expansion of silicon is merely 2.6 ppm/°K and the G and t eff are considered as 
temperature independent, the thermal drift efficient of current related sensitivity can be given 
by [12,13]: 

1 dS, 

asi = ~S~~dT =CCrH ~ aN (19) 

where, a rH and a N are the temperature coefficient of the Hall factor and the carrier concentration of 
N-well, respectively. For a N-well doping of 4 x 10 cm , a rH increases from 0 to 500 ppm/°K in a 
industrial temperature range (240 K-400 K) [12]. With a lower doping of N-well, a rH is almost 
independent of N-well doping and the temperature, showing an approximately 700 ppm/°K constant 
value [12]. On the contrary, a N decreases from 500 ppm/°K to 50 ppm/°K for the N-well doping of 
4 x 10 cm in the same temperature range [13]. As a result, a N could just right compensate a rH at 
the room temperature, and zero temperature coefficient of sensitivity could be obtained. Considering 
the temperature dependency of a rH and a N , the drift of the current related sensitivity is about in the 
range of +500 ppm/°K throughout the industrial temperature range. For the other doping values of 
N-well, a rH and a N can be obtained by referring to relevant data. Thus, considering this thermal drift 
effect of the sensitivity, Equation (20) can be rewritten as: 

S I (T,U pn ) = S I (U pn )[l + a SI -(T-300K)] (20 ) 

2.4. Piezo-Resistance and Piezo-Hall Effects 

When a Hall plate is assembled, its performance is deteriorated by two physical stress-related 
effects, i.e., piezo-resistance and piezo-Hall effects. The piezo-resistance effect due to packaging stress 
provokes relative variations of N-well resistances. In combination with technology variation, junction 
field effects and temperature drift, it is the main source of offset. The sensitivity is also impacted by 
packaging stress. The variation of sensitivity is called piezo-Hall effect. Considering this effect, 
Equation (20) can be rewritten as [14,15]: 

S, (a, T, U pn ) = S I (T, U pn )[1 + P l2 ■ (a x + a y )] (2 1 ) 

where a x and a are the mechanical stress in the plane parallel to the Hall plate surface, P 12 denotes 
the piezo-Hall coefficient tensors in x-y plane, which is estimated at 40 x 10 Pa for the N-well 
(4x 10 16 cm" 3 ) [12]. 

Since the mechanical stress changes with temperature, the temperature coefficient of sensitivity 
illustrated in Equation (19) can be rewritten as: 

a SI = a rH ~ a N ~ a piezo-Hall (22) 

For a plastic packaging, the temperature coefficient related to piezo-Hall effect can be defined by [13]: 

a piezo-Hall = ~ Q ' ^12 ' ^ pg ' ( a pg ~ a silicon ) (23) 
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where E is the modulus of elasticity of molding compound, a pg and a silicon is the thermal 
expansion coefficients of molding-compound and silicon, C x is a designed-dependent geometric 
constant. For a typical plastic package such as TSSOP, we can get C l ~ 6 [13], so the approximate 
a P iezo-Haii value of -450 ppm/°K can be estimated [12]. 

3. Comparing Results of Simulation with Experimental Results 

The new simulation model code has been written in behavioral Verilog-A language and was tested 
on a Cadence Spectre simulator tool using AMS 0.8 urn CMOS technological parameters (shown in 
Table 1) [12]. The finger width and finger length of the cross-shaped Hall plate are designed to 40 jam 
and 60 jam, respectively. 

Table 1. Model parameters (using AMS 0.8 jam CMOS process) [12]. 

Parameters Definition Default value 
N d ,nw Doping in N-well 4 x 10 16 cnT 3 
N A ,p+ Doping in top P+ layer 1 x 10 20 cnT 3 
Na,sub Doping in substrate 1 x 10 16 cm 3 
tttw Depth of N-well 4 urn 
t P+ Depth of top P+ layer 0.3 um 
fi n Electrons mobility 950 cm 2 /V.s 
jx u Hall mobility 1,200 cm 2 /V.s 

Rtci First temperature coefficient of N-well resistance 1%/°K 

Rjci Second temperature coefficient of N-well resistance 20 ppm /°K 

P ]2 Piezo-Hall coefficient 40 x 10 11 Pa 1 

asi Temperature coefficient of S, +500 ppm/°K 

oipiezo-Haii Temperature coefficient of piezo-Hall effect -450 ppm/°K 

To show the correctness and accuracy of this model, the corresponding experimental results of the 
Hall plate fabricated using the same technology given in the literature [12] are compared with the 
model's simulation results. First, we performed the simulation of Hall voltage versus magnetic field at 
room temperature. The simulated and test results when the input bias current is 1 mA and the Hall 
plate is liberated of packaging stress, are plotted in Figure 5. It can be observed that the simulated Hall 
voltage is proportional to magnetic field intensity. When the magnetic field intensity is 2.5, 5, 7.5, 10, 
12.5 and 15 mT, the simulated Hall signal is 0.172, 0.345, 0.517, 0.69, 0.862 and 1.035 mV, 
respectively. While the value of the measured Hall voltage is 0.185, 0.37, 0.551, 0.735, 0.896 and 
1.058 mV for the corresponding magnetic field intensity. We can see that the simulated current-related 
sensitivity is 69 V/AT, while the tested result is 75 V/AT. A very good agreement is thus obvious in 
Figure 5. If a mechanical stress in the CMOS Hall plate is estimated at a x = a y = -70 MPa for a typical 
plastic packaging [14], which will lead to a variation of the simulated magnetic sensitivity of about 
5% at room temperature compared with the stress-free sensitivity. Secondly, the simulation of the 
N-well sheet resistance versus input voltage was implemented at room temperature without 
packaging stress influence. 
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Figure 5. Comparisons between the measurements and the model simulation for the output 
Hall voltage with 1 mA biasing. 
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Figure 6 illustrates the comparison between the simulated and tested sheet resistances per square 
versus variation of input voltage (sweeps from 0 V to 5 V). The measured sheet resistance per 
square changes from 493 QJy to 648 H/y. By comparison, the simulated sheet resistance per square 
changes in the range (506 QJy-622 QJy) with a small error. In addition, the characteristics of 
magnetic sensitivity vs. temperature drift were simulated. 

Figure 6. Comparisons between the measurements and the model simulation for the sheet 
resistance of N-well dependence of input voltage. 
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The measured and simulated relative variations of the current-related sensitivity related to the 
value at room temperature as a function of temperature for the zero-stress mounting of the Hall plate 
is demonstrated in Figure 7. In this temperature dependence of sensitivity simulation, we assume the 
zero temperature coefficient a SI of Hall plate takes place at 27 °C, and a SI linearly changes from 
-500 ppm/°K to +500 ppm/°K in the temperature range from -40 °C to 1 10 °C. It is obvious that a 
good accordance is achieved between simulation and experimental results for a die absence of 
packaging stress influence. Meanwhile, both simulated and tested results of the thermal drift of S t 
influenced by the plastic packaging stress as a function of temperature are also shown in Figure 7. 
The simulated thermal drift of S, for a typical plastic package (TSSOP) is also in good agreement 
with the measured results. It should be pointed out that the piezo-resistance and piezo-Hall effects in 
packaged Hall sensors are very complex issues, which cannot accurately be modeled by only a small 
number of key physical and technological parameters. Therefore, a very accurate simulation result 
cannot be achieved in some cases. 

Figure 7. Comparisons between the measurements and the model simulation for the 
relative variation of the current-related sensitivity as a function of temperature. 
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Finally, the ac simulation of the Hall plate was performed at 3 V DC bias. The ratio of finger 
length to finger width is fixed to 1, while the finger length is taken as a parameter, changing from 
40 um to 120 um with a step of 40 um. The simulation results in Figure 8 show that the smaller Hall 
plate has the higher corner frequency and the -3 dB bandwidth highly exceeds the one MHz range 
for the largest Hall plate, indicating that any limited frequency response of Hall plate within the 
working range (usually below 1 MHz) cannot be observed. 
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Figure 8. Ac simulation gain versus frequency for three Hall plates with different dimensions. 
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4. Conclusions 

An equivalent circuit simulation model for a CMOS-integrated Hall plate has been improved. The 
structure of the model consists of a passive network, including eight non-linear resistances, four 
depletion capacitances and four current-controlled voltage sources. The model completely takes into 
account the non-linear conductivity effects, geometrical effects and temperature effects. Meanwhile, 
the packaging stress influence on Hall plates is also considered to a certain degree. In addition, the 
model only needs a small number of key physical and technological parameters. The model has been 
implemented in Verilog-A hardware description language and was successfully tested with the 
standard EDA tool Cadence. For testing the model correctness and accuracy, the model simulation of a 
Hall plate were performed using AMS 0.8 um CMOS technology parameters and are compared with 
the measured results reported in the literature [12]. A very good agreement is obtained. It should be 
noted that if those key technological parameters such as N-well sheet resistance, Hall mobility, etc. can 
be calibrated by the measurements of the Hall sensors fabricated in a standard CMOS technology line, 
more accurate model simulation results could be achieved. 
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